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NATIONAL ADVISORY COMMITTEE FOE AERONAUTICS 
RESEARCH MEMORANDUM 

PRELIMINARY CORRELATION OF EFFICIENCY OF AIECEAFT GAS-TURBINE 
COMBUSTOES FOE DIFFERENT OEEEATING CONDITIONS 
By J. Howard Childs 


SUMMARY 

An empirical correlation, is presented of the combustion- 
efficiency data previously obtained in experimental investigations 
with 14 aircraft gas-turbine combustors. The precision of the cor- 
relations is not good for some of the combustors. The correlations 
are of value, however. In predicting the efficiency at different 
operating conditions and in comparing the performance of different 
combustors from data obtained in unrelated tests. 

A theoretical analysis based on the kinetics of a bimolecular 
chemical reaction is also presented which, with additional assump- 
tions regarding the mechanism of the combustion process, yields the 
correlation parameter previously obtained by empirical means. The 
theoretical analysis yields equations that predict the shape of the 
curves obtained by plotting the experimental data. 


INTRODUCTION 

Investigations of the performance of various gas-turbine com- 
bustors at the NACA Lewis laboratory have yielded a large amount of 
combustor performance data, but no method has been presented for 
correlating these data in terms of fundamental parameters. Refer- 
ence 1 showed that combustion efficiency and maximum obtainable 
temperature rise were decreased by a decrease in combustor- inlet 
pressure or temperature and by an increase in combustor- inlet 
velocity. No parameter was expressed, however, for combining these 
effects to make possible a prediction of combustor performance for 
any given operating conditions. 

A parameter was derived by empirical methods which served to 
correlate these data, and the correlations are presented herein. 

A theoretical analysis of the combustion process was also made to 
illustrate - that the empirically-derived parameter could also be 
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derived from a theoretical approach. This work was carried out 
simultaneously with the experimental phases of the research. The 
theoretical analysis and methods of correlating combustor perform- 
ance data that are presented in this report were arrived at as 
early as 1944. Attempts to refine the analysis and the use of 
entirely different theoretical analyses have failed to yield a 
better means of correlating the experimental data. The method of 
correlating combustor performance data that is presented herein 
is therefore as good as could readily be achieved with presently 
available experimental data and with the present knowledge of the 
combustion process. 


ANALYSIS 

Empirical. - The data of reference 1 show that the combustion 
efficiency is affected by combustor- inlet pressure, temperature, 
and velocity in a maimer that later investigations proved to be 
characteristic of all gas-turbine combustors. In reference 1 com- 
bustion efficiencies are recorded when each of these three variables 
is independently altered. With these data, the log of the combustion 
efficiency is plotted against the log of each of the three variables. 
The slopes of the resulting lines show the manner of variation of 
the efficiency with each of these variables . The significant para- 
meter that resulted from this empirical analysis is T?^S i /v r . In 

this parameter, P^ is the combustor- inlet static pressure in pounds 
per square foot absolute; T i , the combustor- inlet static tempera- 
ture in °E; and V r , the combustor reference velocity in feet per 
second computed from the total air mass flow rate, P^, T.^, and 

the maximum cross-sectional area of the combustor flow passage. Com- 
bustion efficiency is correlated as a function of PjT^/Vj. . 

Theoretical. - A further analysis was made to illustrate that 
the data correlations could be predicted from theoretical consider- 
ations . 

A theoretical treatment of gas-turbine combustion is difficult 
because of the many different processes that occur in the combustors. 
The fuel must be vaporized, mixed with air, ignited, and oxidized to 
the final products of combustion. These combustion products must 
then be mixed with dilution air to reduce the temperatures to values 
that can be tolerated by the turbine blades. The combustion can be 
visualized as a competition between the., conversion processes 
(vaporization, mixing, ignition, and oxidation) and the quenching 
that occurs when the reacting mixture is cooled and diluted by the 
dilution air and when the reacting mixture comes in contact with the 
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relatively cool walls of the combustor liner. If the rate of any one 
of the conversion processes is substantially less than the rates of 
the others, this one process will govern the over- all rate an d hence 
vill determine the combustion efficiency. The relative rates of the 
conversion processes can be espected to change as combustor operating 
conditions are changed, and thus it may be that entirely different 
processes or different combinations of processes will deter mine the 
combustion efficiency at widely different operating conditions. 


The theoretical analysis was based on the assumption that the 
chemical reaction (the oxidation of the fuel) constitutes an over- 
all rate-determining step in the combustion. This oxidation probably 
occurs by a chain mechanism. It is frequently true, however, that 
the kinetics of chain reactions are dependent entirely upon the 
kinetics of a single, slowly occurring reaction within the chain. 

The analysis was therefore based on the kinetics of a bimolecular 
chemical reaction. Additional simplifying assumptions were required 
in order to obtain from this analysis the parameter pre- 

viously derived by empirical methods. The theoretical ’analysis is 
presented in detail in the appendix, where the symbols used herein are 
defined. The final equations resulting from that analysis are as 
follows : 
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where and Ng are the concentrations of the two reactants in the 

bimolecular chemical reaction; and Eg, Kg, K-j_q, and K-^ are 

constants. Equation (l) applies for a given combustor operating with 
a single fuel at a fixed fuel-air ratio for the general case where 
N^< N-g. Equation (2) should apply for the same conditions of oper- 
ation only if the combustion occurs in a zone where = Kg. 

Equation (3) is a general statement of either of the preceding two 
equations . 
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The arbitrary assumptions on which the theoretical analysis is 
based are as follows: 

(1) Some bimolecular chemical reaction is the over-all rate- 
controlling step in combustion. 

(2) The temperature T in the reaction zone is constant. This 

assumption renders constant the important kinetics term e"* 1 / 51 , 
where 1 is the energy of activation for the chemical reaction and 
E is the gas constant. 

(3) The volume of the reaction zone does not change with a 
change in operating conditions. 


BESULTS AMD DISCUSSION 

Figures 1 to 14 show values of combustion efficiency tj for 

14 different aircraft gas-turbine combustors plotted against the 
parameter accordance- with equation (3). The parameter 

P^T^/V is equivalent to a dimensional constant times the parameter 

p 

/W a where W a is the air mass flow rate. These combustqr perform- 
ance data were obtained from MCA experimental investigations with 
these combustors. The different combustors are given letter desig- 
nations in figures 1 to 14, and most of the combustors are not 
described herein. Several combustors of both the annular and tubu- 
lar types are included, and some of these are from current production 
engines. The data for figure 1 were taken from reference 1 by reading 
values of from the faired curves (fig. 8 of reference 1) at each 

of four fuel-air ratios (0.011, 0.013, 0.015, and 0.017). The data 
for the other combustors (figs. 2 to 14) were computed from the ori- 
ginal test data, using each individual data point. The data of 
figures 1 to 14 are therefore for a range of fuel-air ratios. 

The precision of the data correlations for most of the combustors 
is not good. Some scatter of the data points is to be expected, 
however, inasmuch as it is seldom possible to accurately reproduce 
values of in day-to-day operation at the same test conditions. 

In repeated check tests at identical operating conditions, values of 
differing by 4 percent are encountered with most combustors, 
and at very severe operating conditions the values of have 

been observed to differ as much as 10 percent with many combustors. 

For the combustor of reference 1, values of may differ as much 

as 17 percent at identical operating conditions owing to the 
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tendency of this combustor to follow two different operating curves . 
.(See figs. 8(b) and 8(c) of reference 1.) It has not been estab- 
lished that other combustors follow such dual operating curves, but 
it is suspected that some of them do. A correlation of experimental 
values of tj^ that for most combustors gives a max imum scatter of 
±2 percent at favorable operating conditions and ±5 percent at 
severe operating conditions is therefore as good as can possibly be 
achieved. A criterion for a good correlation might be ±3 -percent 
maximum data scatter at favorable operating conditions and ±7 -percent 
at severe conditions. In view of these considerations, the data 
correlations of figures 3 and 4 appear good, inasmuch as only a few 
data points exceed the allowable ±7 percent deviation at severe 
operating conditions. The data scatter in figures 7 and 8 is so 
great, however, that the correlation is of little value. 

When the parameter P i T i /V r is used to correlate combustor 
performance data, a separate curve might be expected for each fuel- 
air ratio. Most combustors, however, give a substantially constant 
efficiency for a range of fuel-air ratios (this is evident in the 
curves of fig. 8 of reference 1) and thus the experimental data for 
most fuel-air ratios fall on a single mean curve in figures 1 to 14. 
At very high and very low fuel-air ratios, the efficiencies are 
generally lower than at intermediate fuel-air ratios. The parameter 
P’i^i/'^r does not correct for this effect, and' some of the data 

scatter In the figures is due to this phenomenon. 

For each of the combustors the curves of figures 1 to 14 have 
the same general shape. At high values of the parameter P^T^/Vj., 

combustor performance is quite satisfactory; the combustion effi- 
ciency Is high and is not very sensitive to changes In operating 
conditions, as evidenced by the gradual slope of the curves. At 
low values of 3?iTi/V r , however, combustor performance is not 

satisfactory; the efficiency is low and decreases rapidly as 
JP^i/Vr decreases. This characteristic shape of the performance 
curves leads to the concept of a critical value of P^T^/7 r , 

which separates the satisfactory from the unsatisfactory range of 
operating conditions. For each combustor the critical value of 
■ P I T i/ V r is artli ' fcraril -7 selected as the value corresponding to 

= 0.85 on the curves of figures 1 to 14. 
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These critical values of can use< ^ ^ orm a 

rating scale that iri.ll indicate the relative abilities of dif- 
ferent combustors to perform satisfactorily at severe (lov P-j_Tj_/v r ) 
operating conditions. The better combustors vill have lower values 
of critical P i T i /7 r . Table I presents such a rating scale with the 
14 combustors listed in order of increasing critical ■'■k® 

combustion efficiency of a gas-turbine engine at altitude can be 
improved without an attendant increase in engine frontal area by 
replacing the combustor with another combustor that is higher on 
the rating scale in table I. 


Inspection of figures 1 to 14 shows that for values of P^T^/V r 
below the critical (i^C 0.85), the slopes of the curves of 
against P^T^/V r are in general greater for those combustors having 
lower critical values of P^T^/V r . This trend might be expected; 
absolute pressure limits of inflammability exist, and the combustors 
of better design merely extend the range of P^T^/v r for which they 
give high values of closer to these limits. 


Equations (l) and (2) are of different form; equation (l) holds 
for the general case where E A < Eg, and equation (2) should be valid 

only for the special case where E A = Eg. Figure 15 shows the data 


of figure 3 plotted so as to obtain the straight-line relation 
expressed by equation (l); that is, the log of the function 

/ 1 - n, 1/(1 - T], ) is plotted against the parameter P.T./V 

\ h b 7/ 

for two arbitrary values of the ratio E A /Ng ( N A / N B = V 2 and - lA) • 


The predicted straight-line correlation is obtained with either 
value of the. ratio n a A b * Figure 16 shows these same data plotted 

in accordance with equation (2). The data do not fall on a straight 
line; equation (2.) therefore does not apply and hence the reaction 
does not occur in a zone where N = E , which is in accord with 
probability considerations. 


The fact that the straight-line relation predicted by 
equation (1) is actually obtained in figure 15 means that equation (l) 
is of the proper form and predicts the shape of the experimental 
curves; thiB lends credence to the theoretical analysis that yielded 
equation (1). It is emphasized, nevertheless, that the theoretical 
analysis is not of such a nature as to establish definitely chemical- 
reaction kinetics as the over-all rate- controlling factor in the 
combustion. Several arbitrary assumptions were required in order to 
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arrive at the parameter • It might be possible to derive a 

parameter that will correlate the data by means of a theoretical 
analysis entirely different from that presented herein; such a 
parameter, to be successful, would have to be quite similar to 
P^T^/Vr ) and the equations would have to be a form of the second- 

order law like those derived herein. 

Attempts have been made to correlate combustor performance 
data by using parameters derived from each of the assumptions that: 
(l) fuel vaporization, (2) mixing by molecular diffusion, (3) 
mixing by turbulent eddy diffusion, and (4) quenching by wall 
contact control the over-all rate of combustion. The parameters 
based on each of these assumptions gave no correlation whatsoever 
of the experimental data. Attempts have also been made (using 
empirical methods) to refine the equations presented herein so as 
to obtain better correlation of the data; no substantial improve- 
ments in the correlations have been achieved, however. It is 
therefore concluded that the data correlations presented herein 
are as good as can be readily obtained with the present knowledge 
of the combustion process and with presently available performance 
data. 


Data that establish the importance of fuel atomization in the 
combustion process are available, as, for example, ref erence 2. These 
data indicate that vaporization of the fuel may play a rate- 
determining role. The differences in fuel atomization resulting 
from the use of different injection nozzles may serve primarily to 
alter the volume of the combustion zone, however, and thus such 
data do not necessarily invalidate the analysis presented herein. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, • 
Cleveland, Ohio. 
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APPENDIX 

THEORETICAL ANALYSIS OF TURBINE-ENGINE 
COMBUSTION PROCESS 

Fuel atomization, dispersion, and vaporization, mixing of the 
fuel and primary air, ignition, oxidation of the fuel by a chain 
mechanism, and quenching and dilution by secondary air are important 
phases of the combustion process in gas-turbine combustors. For 
certain ranges of operating conditions it is possible that the oxi- 
dation process, or more specifically one step in the oxidation 
chain, governs the over-all rate at which combustion proceeds. The 
analysis is therefore based on the kinetics of a bimolecular reac- 
tion; that is, a reaction involving the collision of two molecules. 
The range of operating conditions for which combustor performance 
data satisfy the equations will indicate the range for which chem- 
ical kinetics, as applied in this particular analysis, may govern 
the over-all combustion rate. 


Outline of Analysis 

A simplified derivation of the necessary kinetics equation is 
first made. This derivation is similar to those appearing in text- 
books on reaction kinetics (for example, reference 3) . Such sim- 
plified derivations involve assumptions and approximations; refer- 
ence 4 presents a rigorous derivation of these equations. The 
kinetics equations are then applied to the case of a typical gas 
turbine combustor, and the variables in the equations are replaced 
by variabies that are measured in combustor investigations. Finally, 
performance data obtained with aircraft gas-turbine combustors are 
plotted in accordance with the derived equations. 


Symbols 

The symbols used in the analysis are as follows: 

A q average cross-sectional area of reaction zone, sq ft 

maximum cross-sectional area of combustor flow passage, sq ft 
energy of activation for reaction, ft-lb/lb 
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ea 


Kj, etc. 


constants 


dw 

d0 


reaction rate, lb mixture reacting/ ( sec ) (cu ft) 


2 distance from upstream end of reaction zone to plane of 

any point vithin reaction zone, ft 

L length, of reaction zone, ft 

n number of molecules per cubic foot, cu ft - '*" 

N number of molecules per pound original mixture, lb - '*' 

P static pressure, lb/sq ft abs. 

E gas constant, ft-lb/(lb) (°E) 

T absolute static temperature, °E 

u mean molecular velocity, ft/sec 

V q velocity based on density, mss flow, and average cross- 

sectional area of reaction zone, ft/sec 

v r velocity based on combustor- inlet density, total air flow, 

and maximum cross-sectional area of combustor flov 
passage, ft/sec 

X fraction of N reacted at any time 9 or at ary location 

vithin combustor 


Z number of molecular collisions of type involved in reaction, 

-3 -1 

per unit volume per unit time, ft sec 

a probability factor, fraction of collisions involving suffi- 

cient energy for reaction to occur vhich actually result 
in a reaction 


\ 


combustion efficiency 


0 time, sec 
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density, Ib/cu ft 

effective molecular diameter, ft 


X 
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P 

a 


cp fraction of collisions involving sufficient energy for 
reaction to occur 

Subscripts: 

A ■ reactant A in reaction A+B — -C+D 
B reactant B in reaction A+B — -C+D 

i combustor inlet 


Simplified Derivation of Kinetics Equations 
For the bimolecular reaction 

A+B— * C+D 

in which A is the reactant present in smaller quantity. The reac- 
tion rate is given by the equation 


dw Zq> 

de ~ Iff a 


(Al) 


The probability factor a will be assumed to be constant, which is 
conventional. 

A molecular collision of the required type (A+B) occurs when 

the distance between the centers of the molecules is (a.+cO/2. 

A B ' 

The average A molecule will therefore sweep out in unit time a 
cylindrical space of length u and having an effective cross sec- 
tion of (jt/4)(G A -*0-g) 2 . The number of collisions of the average A 
molecule with B molecules per unit time will be (je/ 4) (<J A +a-g ) 2 ^ 

Each of the other A molecules is also colliding with B molecules ; 
the total number of A+P collisions per unit time per unit volume 
is therefore 


Z = f (a +a y 

4 A B 


un 


A n B 
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Certain approximations are involved in this derivation, and 
vhen corrections are made for the distribution of molecular veloci- 
ties and other factors, the equation becomes 


Z = ^ ^A +C B ) un A n B 


(A2) 


By substituting for u, PN A (l-X) for n A , and 


N 

PN b (1 - yt X j for V 

s 




N. 


2 " *3 (a A +C %> &*> (\N b )( 1-X) ! - ^ X p‘ 

B 


When — is substituted for p. 


N. 


Z = 


% < Vb t 1 - 2 ) I 1 - jr x ) ^ 


*B 


T e 3/2 t 3/2 

From the Maxwell-Boltzman distribution lav (references 3 and 4) , 


(A3) 


Cf> = 

Equations (A3) and (A4) are substituted in (Al) to obtain 


dv 

d6 


*3 (c A+ a B ) 2 n b (1-x) (i - ^ s) P 2 e- E / m 


3/2 3/2 

R ' T 


(A4) 


(A5) 


Now 


dX 1 dw _ 
dS = p W 


h (W 2 % d-*) (1 - r x ) F ° 

B 


-e/et 


1/2 Bpl/2 


E ' T 


(A6) 


The equations apply to any given sample of the reacting mixture 
as it passes through the combustor. The location of the sample and 
the value of X for the sample vary with time. At the outlet of 
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the reaction zone , 2 will equal because the particular reac- 
tion under consideration is the one governing the over-all process. 
Thus, 


X = when 2 = L and 0 = L/V (A7) 

and 

X * 0 when 2=0 and 0=0 (A8) 

It will now be assumed that the temperature in the zone where 
the reaction A+B— - C+D occurs is substantially constant and inde- 
pendent of the combustor operating conditions and the combustion 
efficiency. Thus, T is constant. With this simplifying assump- 
tion it is possible to separate variables in equation (A6) and 
integrate 



1 - 


N. 


In 


B 


1 -X 


+ K a = 


^ ( W 2 n b m e " E/EE 

E l /2 t 1/2 


(A9) 


Eearranging terms and substituting from equation (A7), X = t], 
when 0 = L/V & , d 


1 — ru 

1 " -.T ^ 


N- 


ln 


'B 


1 - 


% 


% - 


*3 






L Pe 


-e/et 


V a E ' 


1/2 t 1/2 


(A10) 


A value of zero would be obtained for the constant Kg by 

substituting the limits in equation (A8), X = 0 when 0=0, 
into equation (A10) . This relation is not a true boundary con- 
dition for equation (A10), however, because at some operating 
conditions X = 0 when 0 > 0, Flame extinction or blow-out 
occurs whenever the velocity through the combustor exceeds a 
critical value, which is a function of the pressure, inlet tempera- 
ture, and fuel-air ratio in the combustor (reference l). At 
velocities above this critical value the quantity 0 » L/V a > 0 

and yet X = 0 since burning lias ceased. The constant It is 
therefore not equal to zero. “ 
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Application to Typical Gas-Turbine Combustor 

The foregoing equations are expressed in terms of the values 

of P, T, and V in the reaction zone. In experimental investi- 

a 

gations these variables are not measured; the values of the combustor 
Inlet variables P^ , T^ , and V r , are usually determined however. 
The values of P, T, and Y & must therefore be expressed in terms 
of combustor- inlet conditions. 

The pressure loss through a gas-turbine combustor is small com- 
pared with the inlet pressure, so the pressure P in the reaction 
zone is nearly equal P^. In order to make some allowance for the 

difference between P and F^, the following substitution Is made 

P = Ky P ± (A13) 

Actually JL, is not a constant, but it can be treated as such here 
by using a mean value for K_ because its value never differs greatly 
from unity. 

Between the combustor inlet and the zone where the reaction 
A+B — -C+D occurs, the temperature Of the mixture will markedly 
change owing to the processes of fuel vaporization, chemical reac- 
tion (the reactions preceding the particular one being considered 
in the analysis), and heat transfer to and from the surroundings. 

It has already been assumed that the particular reaction under con- 
sideration occurs at a substantially constant temperature so that 

T = constant 


VmL A > vit* tSfi 


(A14) 
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The value of 7 is given "by the relation 

ct 


T 

T a -*8*9Tri®T 


(A15) 


where Eg is the fraction of the total air mass flow that passes 
through the reaction zone , and Eg is the area ratio A /a . 
Equation (A13) is substituted in (A15) to obtain r & 


V. = 


W Y r 

*1 


E- 


(A16) 


Substituting equations (A13) , (A14), and (A16) in 
equation ■ (AlOj , 


1 " N ^ 

In + Eg = E 3 


<W <W • 


It, 2 P A 


It 


E l/2 I 3/2 




(A17) 

On the right side of equation (A17) the terms are grouped into: 

(l) those dependent on fuel type and fuel -air ratio, ( 2 ) those 
dependent on combustor design, and (3) those dependent only on 
the operating conditions. Substituting (A13) , (A14), and (A16) 
in (A12) , 


\ 

i + E c = E* 

1 - V 6 3 


( vV 


N b e _E/ ^ T IE 7 2 P ± T ± 


R l/2 t 3/2 


% 


(A18) 


Equation (A17) is the general solution, and equation (A18) applies 
only if the reaction A+B — -C+D occurs in a zone where K. = 

For a given combustor operating with a given fuel at a filed fuel- 
air ratio, equation (A17) becomes 

H. 

1 " ST \ P.T, 

1,1 TT^~ * h ■ *M <« 


" i»^- LUlJiUEBTiAL-^ 
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For these same conditions, equation (A18) becomes 


% 

1 - % 


+ Eg 


= E- 


P T 
ii 


’ll Y„ 


( 2 ) 


which is applicable 
in a zone where N. 
follows 


only if the bimolecular reaction happens to occur 
= Hg. Either equation can be expressed as 


% 



(3) 


EEFERENCES 

1. Childs, J. Howard, McCafferty, Richard J. , and Surine, Oakley W. : 

Effect of Combustor-Inlet Conditions on Performance of an 
Annular Turbojet Combustor. MCA Eep. 881, 1947. (Formerly 
MCA TN 1357). 

2. McCafferty, Bichard J.: Effect of Fuels and Fuel-Nozzle Char- 

acteristics on Performance of an Annular Combustor at Simulated 
Altitude Conditions. MCA EM E8C02a, 1948. 

3. Getrnan, Frederick H. , and Daniels, Farrington: Outlines of Theo- 

retical Chemistry. John Wiley & Sons, Inc. (New York), 1937. 

4. Eassel, Louis S.: The Nineties of Homogeneous Gas Beactions. 

Chem. Catalog Co. (New York), 1932. 








iHa.iirA#, 












DCHiCl.M*! 








Llllbtlc 










I 





















i i 

N Af A ^ 


























i i " 


. |Q l I l l ■ - l l I -j I I I L . 1 I 

‘ 28 30 32 34 38 38 40 4&XI0 3 


P|T|/V r 

(b) P,T 1 /Y r >£ 8 xl 03 . 

Floure 3. - Concluded. Experimental data obtained eith coebustor C. 

Feel. AM-F-58. 


sssrt 


' !>'i ti- 


r i 


NACA RM E 5 QFI 5 . 




FI dura 4, 





VUN3CLJ 





20 22 24 26 28 30 32 34 38 38 40 42 44 4B 4&CI0 3 

PlT,/V r 

(b) P,T i /V r >20xl0 3 . 


Figure 4. - Concluded. Exporlaenta! data obtained with coobostor D. Fuel, AH-f-28- 




£S£I 



NACA RM E50F15. 










NACA RM E50FI5. 






Figure 7 



clareter coahustor of NACA design (coahustor G). 
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(b) P i T|/V r >20xl0 3 

Figure 9. - Corel Bded. Experimental data obtained wlt^ combustor I. 
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[b) P|T,/V r >20xl<)3. 

Figure 10. - Concluded. Experimental data obtained with coa 
buster J. Fuel, AH-F-B8. 
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(b) P | Tj /V r > 20x l()3; fuel, AN-F-32. 

Figure II. - Concluded. Experimental data obtained with coabustor K. 
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(a) Na/Nr =■ 1/2. 

Figure IS. - Experiments! data obtained with combustor C plotted in accordance with equation (1). 
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Figure 16. - Experimental data obtained with combustor C plotted 
in accordance with equation (2). 
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